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Abstract 

After an historical introduction reviewing the successes and failures for heavy 



> 
(N 

■ quarkonium spectroscopy of the nonrelativistic quark model including {v/cf' 

(N 



corrections, the discussion is widened to include light quark pair creation; i.e., 
dynamical quarks. We find that the simple extension of the QCD inspired 
potential model which includes dynamical quarks first proposed by Eichten 



Oh! 

1^ ■ et al. P] is remarkably successful in accounting both for the observed narrow 



state heavy quarkonium spectroscopy and experimental data on threshold 
production of charmed and B mesons in e'*"e~ annihilation. We have studied 
various different models of light quark pair creation. In addition to the original 
Cornell model, we find that a variant of it can also account for the data. The 
cross section data are mainly inclusive cross section measurements. The two 
models give distinguishably different predictions for the energy dependence 
of production cross sections for individual channels, and measurement of the 
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energy dependence of production cross sections for individual channels may 
distinguish between these models. 



I. INTRODUCTION 

In this talk I would like to discuss three sorts of QCD inspired phenomenological models 
which describe QQ physics at low energies. [0] The simplest is the nonrelativistic naive quark 
model (NR NQM). By naive quark model (NQM), I refer to treatments in which the problem 
is treated strictly as a two body problem. As I remind you below, the NR NQM initially had 
many successes. However, there were also some failures which were put right when {v/cf' 
corrections were taken into account. I will discuss separately two types of (f/c)^ corrections; 
(i) those applied to the NQM and (ii) extension of the NQM to take into account light quark 
pair creation. This is one of the two main points of my talk, namely that light quark pair 
creation effects are (f /c)^ corrections to the NQM and are of the same order as effects due to 
the spin dependence of the QQ force and kinematic (f/c)^ corrections. All these effects are 
important in obtaining good agreement of quark model predictions with experimental data. 
There are two types of light quark pair creation effects to be discussed here; namely, (1) 
those that occur in heavy quarkonium states below flavor threshold (vacuum polarization 
effects) and (2) those responsible for production of charmed and B mesons at low energies. 
The effects of virtual light quark pairs on the narrow QQ states below flavor threshold 
can be included through modification of the two-body QQ potential. If one confines one's 
attention only to these states, then the phenomenological QQ potential obtained by fits to 
data includes both the two-body potential and the effect on the spectrum of virtual qq pairs. 
Such a treatment, however, does not allow for relating the effects of virtual qq pair creation 
with actual qq pair creation as seen in the productions of charmed and B mesons. It is also 
somewhat inadequate in that it is difficult to obtain the configuration mixing necessary to 
account for certain spectroscopic results without explicit inclusion of four-body QqqQ virtual 
states. 
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In section II is a brief summary of the successes, and failures, of the NR NQM; and, 
in section III, the improvement when (f/c)^ corrections and the spin- dependence of the 
QQ interaction are taken into account. Section IV presents some QCD inspired potential 
models of light quark pair creation. Section V outlines the method used to include these 
dynamical quarks in calculating both the spectroscopy of the narrow QQ states below flavor 
threshold and the production cross sections for charmed and B mesons above flavor thresh- 
old. Of the various possible models we studied in detail, two seem reasonably successful in 
accounting for both the observed spectroscopy of the narrow states and the observed thresh- 
old production cross sections in e+e^ annihilation. These models are the Cornell model p| 
and an extension of it studied by Zambetakis [^. Section VI compares predictions of these 
models with available production cross section data. These are mainly includive cross sec- 
tion data. Both models give remarkably good accounts of these data but neither are precise 
fits to the data. This section also presents exclusive production cross sections; for these the 
models predict different energy dependences. Measurement of exclusive production cross 
sections may indicate which of these models is the more correct one. 

II. SUCCESSES AND FAILURES OF THE NONRELATIVISTIC NAIVE QUARK 

MODEL (NR NQM). 

This section is something of an historical account which is incomplete and included here 
to give background to the assertion that corrections to the NR NQM owing to light quark 
pair creation, and the coupled channel mixing it induces, are of the same order as those due 
to spin-dependent forces and other {v/cY corrections. 

A. Successes 

Here are outlined the early outstanding successes of the NR NQM. Using a simple QCD 
inspired potential which behaves at short distances like a one gluon exchange potential and is 
linearly rising at large distances, the model has as parameters the strength of the Coulombic 

3 



part and the slope of the hnear part of the potential; in addition, the masses of the c and 
b quarks. With these few parameters, the model accounts quite well for a large number of 
data: 

Masses of narrow heavy quarkonium states (spin averaged). 
Leptonic decay rate ratios ree{'4'')/^ee{'^), etc.. 

Allowed radiative transition rates (given by < r >/i); e.g., -0' xj, xj ~^ J/4'i 
etc.. 

Allowed Ml radiative transition rates (given by < 1 >/i) ; e.g., J/ip ^ rjc ■ 

B. Failures 

Among the successes above were some failures. These were: 
ree(J/V') too big. 

El radiative transition rate for ip' xo too big by a factor 2. 
Forbidden Ml transitions such as ip' — > rjc observed. 

■0(3770) observed in e+e~ annihilation with mass too close to the '0' to be un- 
derstandable as a ^3Si cc state. 

In addition this simple model failed to account for the observed fine and hyperfine splittings 
in the spectra. The items above, and the fine structure of the mass spectra, are corrected 
by taking relativisitc corrections into account. 

III. RELATIVISTIC CORRECTIONS. 

There are two types of relativistic corrections: (i) {y/cY corrections to the NQM; i.e., 
{v/cY corrections to the nonrelativistic potential model in which nevertheless the model 
remains a two-body treatment of the problem, and (ii) extension of the model to take into 
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account dynamical quarks; i.e., light quark pair production. Type (i) relativistic corrections 
have been discussed in this Conference explicitly by G. M. Prosperi, Yu-Bing Dong, L. 
Fulcher, F. Schoeberl, H. Sazdjian, Kuan-Ta Chao, and extensively in the literature; see,e.g., 
and references cited therein. Spin-dependent forces are also regarded as type (i) (f/c)^ 
corrections because they arise in this order in the Breit-Fermi reduction of a particle exchange 
diagram. Estimates of the spin- dependence of the QQ force have been obtained other ways as 
well; see, e.g., and the contributions of J. M. Ball and F. Zachariasen to this Conference. 
The effects that I would classify as type (ii) relativistic corrections are those that arise owing 
to virtual light quark pairs. These are of the same magnitude as those of type (i) and arise 
in relativistic quantum field theory in order (v/c)'^. Below is a list of relativistic corrections 
which have been found to correct the 'failures' mentioned above. In most cases several effects 
contribute. 

a. spin-orbit QQ interaction 

b. spin-spin QQ interaction 

c. tensor forces 

d. spin-independent (f/c)^ corrections to the QQ Hamiltonian 

e. direct QQ ^Di - photon coupling 

f. light quark pair creation 

g. coupled channel mixing. 

In Table I. are indicated which of the above items tend to correct the so-called failures. 
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put right 


by 


Fine structure splitting 


a 


El rate ip' Xo 


a and g 




g 


Forbibben Ml rates 


d and g 


V'(3770) explained 


c, e, and g 


(T(e^e~ — >■ charmed mesons) 


f 



Table I. 



I have not included hyperfine splittings in the above discussion. This is because I believe 
that remains an unsolved problem. 

The bulk of this paper is about light quark pair creation. The foregoing was included 
to indicate that inclusion of light quark pair creation improves the the agreement of NQM 
results with experimental data in two ways. Not only does it extend the model to describe 
quarkonium decays to charmed and B mesons, it also improves the agreement with the 
narrow state spectroscopic data. 

Before addressing explicitly light quark pair creation effects, it is appropriate here to 
remind you that relativistic effects such as spin-orbit couplings, as well as pair creation, 
reflect the Lorentz transformation properties of the interaction. For example, in a Breit- 
Fermi Hamiltonian the spin-orbit interaction coming from vector and scalar exchange have 
opposite sign and differ by a factor 3. Correspondingly pairs created from vector and scalar 
exchange are produced in 1 and O"*""*" states, respectively. At present it is not clear that 
the spin-orbit interaction is correctly viewed as an effect found in a Breit-Fermi reduction 
of vector and scalar exchange interactions. However, when this assumption is made and 
the result compared with data the observed fine structure splittings indicate that the long 
range (linearly rising) part of the potential comes from scalar exchange. We included light 
quark pair creation in such a model and found this coupling to decay channels incapacitated 
the model. (See below.) On the other hand, models in which the quarkonium potential 
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as a whole is assumed vector exchange can account well for the (spin-averaged) masses 0; 
'vector' models also account reasonably well for observed threshold charmed and B meson 
production cross section data. Therefore, it may be reasonable to assume there is some 
truth ion these 'vector' models. They are not inconsistent with the analyses of Eichten 
and Feinberg |^ and Gromes [Q] who obtain in a expansion of QCD spin- dependent 

QQ interactions similar to those obtained in the Breit-Fermi reduction of a vector short 
range and scalar long range potential. 



IV. INCLUDING DYNAMICAL QUARKS 

Eichten et al. proposed a quarkonium potential model that includes dynamical quarks 
nearly 20 years ago. [0| Their framework, with various different dynamical assumptions, has 
been used by most authors writing on this subject since then. Straightforward inclusion 
of light quark pair creation in a QQ model requires quantum field theory and immediately 
confronts us with a many body problem to solve. Eichten et al., hereafter referred to as the 
Cornell group, simplified the problem to manageable proportions. Before going into details 
about model predictions for charmed and B meson productions, I would like to briefly 
outline their framework. Its structure incorporates fundamental properties of QCD, namely 
confinement and asymptotic freedom. Consider the QCD vacuum polarization diagram 
shown below. 



Q 




Q 



For loop momenta large compared to constitutent quark masses, this diagram can be treated 
perturbatively and is taken into account using the QCD running coupling constant. However, 
diagrams with small loop momenta should be treated nonperturbatively. For low momenta, 
higher order diagrams with multiple gluon exchanges should be taken into account. These 
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were included by the Cornell group invoking duality. They used a phenomenological field 
theoretic model for qq creation and then took the strong interactions into account by replac- 
ing the QqqQ intermediate states by Qq and qQ mesons (bound by the same potential as 
binds QQ ). This reduces the many body problem to a tractable pair of two-body problems. 

First one solves a potential model for bound QQ , Qq and qQ states. The Cornell group 
used 

(1) 

Then, as in the Wigner-Weisskopf treatment of coupled channel problems in nuclear physics, 
one solves a coupled channel two-body problem that can be represented by the Hamiltonian 



H 



V ^ ^QmQ J 



(2) 



which operates in a space of bound QQ states and bound Qq and qQ states. In the Cornell 
model the QQ states are eigenstates of 

Hqq = ^ + V + 2Mq. (3) 

and the Qq and qQ mesons eigenstates of the corresponding Hamiltonian (with relativistic 
correction for the light quark). In the calculations, however, measured (or to be measured) 
values of the Qq and qQ meson masses are used. The piece denoted by Hq^^q in @j is 
the kinetic energy operator for the two-meson states. The off-diagonal piece h couples the 
QQ states to their two-body decay channels. (Final state interactions are neglected; i.e., 
interactions between the bound Qq and qQ mesons are neglected.) 

To complete the framework described above, we need to put in a mechanism for pair 
creation. In the Cornell model, pairs are created by the same interaction that binds; i.e., by 
diagrams such as 




where the corkscrew hne represents the propagator of an (instantaneous) interaction whose 
Fourier transform is the quarkonium potential V. The quark vertices, correspondingly, were 
taken as nonrelativistic reduction of a vector exchange. One may vary these assumptions 
and consider other dynamical models of quark pair creation. For example, one could assume 
scalar rather than vector exchange. As mentioned above, the fine structure in the mass 
spectra indicate that the long distance part of the potential may be due to scalar exchange. 
These two assumptions give quite different results; a vector interaction creates a pair with 
quantum numbers 1 while a scalar creates a pair with the quantum numbers of the 
vacuum. We can report here calculations based upon various different assumptions about 
the dynamics of qq pair production. They are tabulated in Table II. 

h 



Model 



ex 



Cornell 



Q 
— >- 




V 



Zambetakis 



Q 
— >- 




V(r) 



V 



Grotch, Zamb, 
Byers 



Grotch & 
Byers 





-K/a + fC 
+ (1-f)C + r/a' 
{o«f <1} 



order 
QCD 



V 
S 



V 



Le Yaouanc 
et al., ... 



Q 
— >- 



q 'qpc 



Table II. 



The dynamical assumptions represented in Table II. are as follows. First, the Cornell 
model: the original Cornell model neglected pair creation from the Coulombic part of the 
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potential and took pair creation only from the long distance part. This is indicated in the 
third column. In the fourth column is indicated whether the exchange is vector or scalar. 
The second row represents Zambetakis' extension of the original Cornell model to include 
the Coulombic part of the potential in the pair creation matrix elements. The third 
row represents our calculations with Howard Grotch in which we took the Coulombic part 
of the potential to be vector exchange and the long range part to be scalar. There is an 
ambiguity in implementing this idea; namely where does short become long? Quantitatively 
this translates into the question of how to treat the constant term C in ([^). For vector 
coupling, C plays no role; the pair creation matrix elements only depend on the gradient of 
the potential. However with scalar coupling the potential behaves like a mass term and the 
pair creation matrix elements depend on V rather than VV^. Thus the constant C can be 
significant. Fits to the data indicate C is large and negative; of the order of —800 MeV. 
To take this ambiguiuty into account, we introduced another parameter / which specifies 
which fraction of C couples vectorially, and allowed it to vary between and 1. The fourth 
row refers to calculations with Grotch of the contribution of the triple gluon coupling graph. 
We evaluated this graph in Coulomb gauge, and found its contributions were not large; they 
are indeed of order a^. We also studied the other QCD graphs of this order and found 
none seems to be singular. We therefore conclude that such effects are not at this stage 



significant. In the last row is represented the QPC model of Le Yaouanc et al.. |jT2| It is 
a model sometimes referred to as a flux tube or string breaking model. In it the qq are 
produced with the quantum numbers of the vacuum, and the pair creation matrix element 
is a constant which can be included in calculations like the above by taking it to be given 
by a constant potential rriq'yQPc / This model gives interesting results which have been 



discussed in a number of papers; see Ref. [12 



In order to get the matrix elements of h from the diagrams in Table II., one must sandwich 
them between QQ wave functions and Qq and qQ . (Generally the QQ wave functions are 
taken to be eigenstates of (0), and for simplicity the Qq and qQ wave functions are taken to 
be Gaussian fits to bound state wave functions bound by the r/a^ part of the potential.) 
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We report here the results of our studies of first four of these models. Mainly what we 
have to report are our results for the Cornell and Zambetakis models because (1) we found 
the order graphs were not likely to be large enough to be important; (2) the hybrid 
model in which the long range part of the potential was considered scalar exchange was 
unsuccessfull (explained below); and (3) both the Cornell and Zambetakis models were able 
to account for available data - a large amount of data on the narrow states below threshold 
and on production cross sections for charmed and B mesons above threshold. This is will 
be discussed in detail in the last section of this paper. 

We end this section with a brief explanation of why we think that the hybrid model of 
short distance vector and long distance scalar exchange is not a viable model. Allowing / to 
vary between and 1, we could not find values of the parameters in the potential (|I|) that 
allowed a fit to the narrow state spectroscopic data. The reason for this is that in the region 
of r where the quarkonium wave functions are large V changes sign, going from Coulombic 
at short distances to linearly rising at large distances, and consequently, when it is coupled 
as a scalar, its contributions to matrix elements of h are relatively small differences of two 
large numbers. Owing to this the matrix elements of h vary erratically as one goes from 
one of the low-lying quarkonium states to another; consequently their contributions to level 
shifts are unstable against small variations of the parameters. 

V. MODIFICATION OF NQM 

Before discussing our numerical results, it may be of interest to outline here how the 
inclusion of dynamical quarks changes the NQM. This can best be seen in a Hamiltonian 
framework in which the light quarks have been integrated out and one has remaining an 
effective Hamiltonian in the QQ space. This is a peculiar Hamiltonian because (i) it is 
energy-dependent; and (ii) for energies above threshold for meson production, it develops 
an anti-hermitian part. It is, however, quite interesting because it gives mass shifts and 
configuration mixings and also gives production cross sections for charmed and B mesons in 
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e"'"e~ annihilation. 

In the space spanned by eigenstates of (|), hereafter called NQM states, the effective 
Hamiltonian is (see, e.g., Refs. 

HeffiW) = M^^^^ + Vl{W) . (4) 

where M.^"-"^^ is a diagonal matrix whose elements are the eigenvalues of (|^). The Vt matrix 
is second order in h] the explicit expression for it is given below. It is perhaps useful to state 
some more of the properties of i^e// here. 

(1) . The physical masses M^v of the narrow quarkonium states are given by the 
solutions to 

ife//(M;v)a^= M^va"^ (5) 

where denotes a state vector in the space spanned by NQM states. A physical 
state ^'n may be expanded as 

*N = E«fV^^ + E&n'^n- (6) 
i n 

where tpi are NQM statesand (pn are continuum two-meson states. The af are 
the components of the 'eigenvectors' in (^. The problem is not exactly an 
eigenvalue problem. The dependence on Mjv is highly nonlinear. Nevertheless 
solutions can be found. |^ 

(2) . H^ffiW) is not hermitian when W is above the threshold of the lowest mass 
two-meson continuum. The quarkonium states whose masses are above threshold 
are poles in the complex W plane. They are poles of the effective propagator 
Q{W), where 

g{W) ={W - H^ffiW))-' . (7) 
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For W above threshold, Q has the structure 

n{W) = R{W) - iT{W)/2 (8) 

with R and F real symmetric matrices. Unitarity requires the diagonal elements 
of F to be positive. 



The matrix elements of are given by are given by 

The sum over n denotes sum over channels and integration over channel phase space; En 
is the channel energy. The channels are labeled by / which specifies channel spins, flavor, 
masses, etc.. Because of this sum over channels, Q is a sum of matrices 

n{W) ^Y.n^f\W) . (10) 

/ 

Each matrix element of O*^'^^ is an integral which becomes singular when W > mi + m2 
where mi and m2 are the masses of the mesons in channel /. These singular integrals are 
evaluated in the usual way taking W in the upper half plane and then the boundary value 
on the real axis from above. 

For W above flavor threshold, charmed and B meson production cross sections in 
e^e~ annihilation may be calculated as follows. The cross section ratio 

Aii' = I'fr , , (11) 

(7(e+ e — > /x+ j 

in one photon approximation, is given by the dispersive part of QiW). With a high speed 
computer and -f^e// in matrix form (see below), Q is easily calculated by matrix inversion. 
There is no coupling between the cc and bb subspaces, so one deals separately with these two 
subspaces and the calculations for the production cross sections of charmed and B mesons 
are done independently. The charm contribution to the cross section ratio Ai?c can be 
expressed as a trace in the cc subspace; viz., 
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AH. = -'^elT.(W^^U (12) 

where Cc = 2/3, and the matrix W is bihnear in the cc -photon couphng. Nonrelativisti- 
cally, matrix elements of W are given by cc NQM wave functions; viz., 

Wij = (0)^,(0) . (13) 

(Relativistic corrections are important here, particularly for charm; see Ref. [0.) In their 
original papers, the Cornell group showed that 

g - =g\Q-Q^)g. (14) 

Thus, from (|12|), (p!^), and (^^, one sees that the exclusive production cross section for 
channel / is 

727r 

Ai?(^) = -— el Tr( W g^ M^f^ g U ■ (15) 

Similarly, one calculates production cross sections for B mesons. From the above it is clear 
that this is a nonperturbative treatment of dynamical quarks. 

With high speed computers it is straightforward to solve, as outlined above, for both the 
narrow states below threshold and the (resonant) cross sections above threshold. One first 
solves d^) for the NQM wave functions and masses. The matrix elements of the VL matrix 
and (^ can then be explicitly evaluated. These matrix elements are discrete because, owing 
to confinement, the entire QQ state space is spanned by a discrete set of eigenf unctions. 
With these matrices one can calculate all of the above quantities. Most computers are now 
able to work with matrices of (almost) arbitrarily high dimension. However, since we are 
interested only in the low lying quarkonium states, we do not need to calculate matrices of 
high dimension. The highly excited NQM states make negligible contributions. 

On the other hand, another approximation involved in a numerical evaluation of fi(W) 
may cause significant error in predicted masses of quarkonium resonances in the continuum. 
In calculating fi(W) one generally approximates the infinite sum over channels / by a fi- 
nite sum neglecting all channels with thresholds greater than some minimum energy Emin 
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which is greater than the maximum value of W for which one is calculating the f2 matrix. 
Though the contribution from any one of these neglected channels is small, their cumulative 
effect can be appreciable because all these channels contribute coherently (negatively) to 
the diagonal elements of Re Q; c.f., Eq. (Cancellations are likely in contributions to the 
off-diagonal elements owing to random phases.) These cumulative mass shifts are small for 
W substantially below Emin', however, they may become appreciable when W comes close 

to Emin- Q 

VI. PRODUCTION CROSS SECTION CALCULATIONS AND COMPARISON 

WITH EXPERIMENTAL DATA 

In this section we report our results on threshold production cross sections, inclusive and 
exclusive, for charmed and B mesons in e'^e~ annihilation and compare them with available 
data. There are two sets of results; one obtained from the original Cornell model and the 
other from Zambetakis' extension in which pair creation occurs at short distances also. The 
available data are mainly inclusive cross section measurements. The agreement with data 
is about the same for both although it differs in detail in the two cases. The two models 
give significantly different results for exclusive channel cross sections. Measurement of these 
may distinguish between them. The degree of agreement of these model predictions with 
data is, in our view, significant because the calculated cross sections are obtained without 
adjustment of parameters or introduction of additional parameters. The parameters are 
the parameters of the potential k, C, and We calculated the masses of the narrow 

^In our calculations we included only the channels with ground state pseudoscalar and vector 
mesons. This is a good approximation for the narrow states below flavor threshold; however, for 
values of W in the continuum we found it necessary, as did the Cornell group, to compensate for 
the neglect of excited meson decay channels by putting in a phenomenological mass shift to fit the 
observed resonance energies • 
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charmonium states using and determined these three parameters by fitting the observed 
masses of J/ip, Xcogi and tp' . Parameter values which fit these data ar^ 



model 


l/a2 in GeV^ 


C in GeV 


k(cc) 


Zambetakis 


0.31 


-0.97 


0.49 


Cornell 


0.22 


-0.85 


0.52 



Using these parameters, we calculated (spin averaged) T masses, charmonium and T leptonic 
widths, etc. and found good agreement with measured values Then we calculated 

charmed and B meson production cross sections above flavor threshold evaluating (|l^) with 
the same parameters. We calculated charmed meson productionfrom threshold to 4.5 GeV, 
and for B from threshold to 11.1 GeV. In the numerical calculations we truncated the cc state 
space after the 4S and 2D NQM states, as in the original Cornell calculations. We found 
that including more states did not signiflcantly alter the results. For B meson production, 
the T(4S) is the flrst hb resonance above threshold; we truncated after the 7S and 4D NQM 
states. Results were stable against increaseing the number of such states. However, we 
also found that the results changed little when we omitted the D-states and only changed 
noticeably near 11 GeV when we omitted the 7S state. The results presented here were 



^ After a detailed and extensive analysis, the Cornell group used the following constitutent quark 
masses Mc = 1.84 GeV, = 5.17 GeV, and for light quarks, m„ = rrid = 335 MeV, and rus = 
450 MeV. We used the same values in our calculations. 

^The slope of the linear potential l/a^ is significantly greater in Zambetakis' model because the 
short distance pair creation increases the coupling of cc to inelastic channels This stronger coupling 
to inelastic channels increases the separation of charmonium energy levels, and requires increased 
strength of the confining potential to bring the mass differences back to their observed values. 

^For the T masses and widths, and for B meson production cross sections, the QCD logarithmic 
decrease of ag was taken into account by using for the values 0.46 and 0.48 for the Zambetakis 
and Cornell models, respectively. 
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obtained truncating bb states at the 6S. The reason it was not necessary to include D-states 
is that owing to the heavier b quark mass, the S-D mixing and photon-D state couphng are 
smaller than in the charm case. 



A. Charmed meson production. 



To understand the complicated energy dependence of the inclusive charm production 
cross section above threshold (3.73 GeV), it is useful to consider the singularities of the 1 
cc propagator Qcc{W) in the complex energy plane; Q{W) is given by (0). The singularities 
are cuts along the real axis and poles in the lower half plane. There is a cut for each open 
channel with a branch point at its threshold. The poles correspond to resonances in the 
cross section. The branch points and poles are shown in Fig. 1 where the positions of the 
singularities are determined by the experimental values of the masses and widths of the 



indicated states the distance below the real axis of each pole is the measured total 
half-width. The separations of the branch points and positions of the poles are to scale. 

D*D* Do-Df+DfDo" D0-D0-+ Dq-Do 




'so 



V(1D) 



V(3S) 



V(4S) 



V(2D) 



Figure 1. 



The corresponding cross section measurements and theoretical curves for R are shown in Fig. 
2. The first resonance above threshold is understood, from the point of view of the Cornell 
model, as the ipllD) though it is not a pure state but has admixtures of other nearby 1 
S and D states. Though this is a relatively narrow resonance as indicated by its proximity 
to the real axis in Fig. 1, it is not as prominent as the S-state peaks, '?/'(4040) and ip{4415), 
owing to the fact that nonrelativistically D states do not couple directly to the photon. As 
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it is primarily a D-state, the '?/'(3770) owes its presence in an e~^e~ annihilation cross section 
to the fact that it has some configuration mixing with nearby S states ( particularly the 2S) 
and a direct coupling to the photon in order a (f/c)^. Similarly the iIj{2D) appears in the 
cross section as a shoulder on the ip{SS) or ?/^(4040) peak owing to its admixture of 3S-state 
and direct coupling to photon. 
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Figure 2. 



The data points in Fig. 2 are from the Review of Particle Properties. [|T^| The curves 
are calculated using the Zambetakis and Cornell models with the contribution to R from 
other than charmed hadrons taken from the data to be 2.65; i.e., the curves in Fig. 2 are 
ARc{W) + 2.65.0 What is remarkable about these models is the degree to which the energy 
dependence and overall normalization of the calculated curves are in agreement with the 
data. Neither model gives a perfect fit to the data; but considering the simplicity of the 
models, and the complexity of the physics, it seems remarkable that the models fit the data 
as well as they do. The Cornell model seems to fit the data in the region of the '?/'(4040) 
better and the Zambetakis model seems to do better in the region of the '?/'(4415). It is 
difficult on the basis of these inclusive cross section data to say if one or the other of these 
models is closer to being correct. 



^ A phenomenological mass shift matrix was used as compensation for the neglect of excited meson 
decay channels; we only included the ground state psuedoscalar and vector meson decay channels 
in our calculations, c.f., section V. 



The models are more easily distinguished in their predictions for the exclusive cross 
sections. In Figs. 3 and 4 are shown the individual channel contributions to Ai?c. 
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Figure 3. 

The full curve is the contribution to Rc from DD channels, both charged and neutral; the 
dot-dashed curve from DD* + D*D; and the dashed curve from D*D* . 
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Figure 4. 



The cross sections for Dg and production are roughly 1/3 of those for non-strange D 
mesons. These are shown in Fig. 5. 
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Figure 5. 

In Fig. 5 tlie full curves are Cornell model predictions for DgDg, DgDl + D*Ds, D*dI., and 
the dashed curves Zambetakis model predictions. The curves for the different channels can 
be distinguished by the fact that they start at different thresholds. Note the scale change. 
This suppression of the strange mesons is natural to these models. It comes from the rris 
being greater than tjIu- The suppression is due, essentially, to two effects: (i) the pair 
creation amplitude is inversely proportional to the light quark mass, and (ii) the strange 
quark meson states are smaller (more tightly bound). 

B. B meson production cross sections. 

The singularity structure of Q{W)g^ relevant for our purpose is shown in Fig. 6. It 
appears simpler than Fig. 1 because the B* — B mass difference is smaller than D* — D and 
because we have indicated only the S-state resonances. In the T system both S-D mixing 
and direct photon- D-state coupling are small. These are effects of order {v/cY and therefore 
smaller than in charmonium. Consequently D-state contributions are not significant in 
e^e~ annihilation cross sections. 

BB*+ B*B BsB| + B^Bs BB** + B**B 
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Figure 6. 

Corresponding to this in Fig. 7, we show the inclusive B meson production cross section 
data |]l4l along with the calculated curves. 
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Figure 7. 



Radiative and beam spreading corrections have been applied to the theoretical curves in 
Fig. 7 so that they can be compared with the data which are uncorrected for radiative 



Owing to the fact that the T(4S) can only decay to BB, it is a relatively narrow resonance 
and prominent in the data. The T(5S), on the other hand, can decay to all nine ground 
state psuedoscalar and vector mesons and is relatively broad. In Fig. 6 its position is far 
from the real axis and relative to its width, its mass is near many sharp thresholds; the mass 
and width are 10.865 GeV and 110 MeV, respectively. Because of these complications, the 
T(5S) does not have a normal Breit-Wigner shape and is hard to see in the cross sections. 
The T(6S) mass is quoted in Ref. [|l^ as 11.019 GeV. We have included in Fig. 7 data only 
out to 10.99 GeV because we think that our calculated values in the region of the 6S are 
unreliable owing to the fact that we have not included channels with B mesons in excited 
states - either with £ 7^ or excited i = mesons; and these are likely to have thresholds in 
the vicinity of the T{QS). 

The curves in Fig. 7 have been calculated without any adjustment of the parameters of 
the model. It is, therefore, in our opinion remarkable that they agree with the data as well 
as they do. As indicated at the end of section V, a mass shift matrix has been used to fix the 



corrections and beam spread.^ 



^We thank Dave Besson, Persis Drell and Elliot Chu for applying the radiative and beam spread 
corrections to the calculated cross sections. 
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peak values of the cross section at 10.580 GeV and 10.865 GeV. Aside from this, the energy 
dependence and the overall normalization of the calculated curves are determined by the 
models. The Zambetakis model (full curve) appears to fit the data around the T(4S) better 
than the Cornell model (dashed curve). However, in the region of the T(5S) the Cornell 
model may be a better fit. 

As in the charm case, the energy dependence of the exclusive cross sections are quite 
different for the two models. The contributions to the cross section ratio R from the various 
channels are shown in Figs. 8 and 9; the full curves are for the BB channels, neutral plus 
charged; the dot-dashed curves are for BB + BB*; and the dashed curves are for B*B . 
Notice that at the higher energies where these curves tend to flatten out they go over toward 
the proportion 1:4:7 which are the spin ratios, discussed in Ref. 0, for BB , BB* + BB*, 
and B*B* , respectively. In Figs. 8 and 9 we omitted the T(4S) resonance region because 
the cross section there is so much bigger than it is in the rest of the range. 
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Finally in Fig. 10 we show the contributions to R from Bs channels; the full curves are 
Cornell model predictions and dashed curves Zambetakis model. The curves for the various 
channels can be distinguished by their thresholds. 
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Figure 10. 

One interesting feature of the curves in Fig. 10 is that both models predict relatively copious 
B,B^ + B,B* production in the T(5S) resonance region. (These predictions are sensitive to 
the Bg mass and the B* — Bg mass difference.) As in the charm case, however, these models 
predict that the strange mesons arc produced at a level of about 1/3 that of the nonstrange 
mesons. Note the change of scale in Fig. 10. 



VII. CONCLUSIONS 

The fact that both the the Cornell and Zambetakis models, as simplistic as they are, can 
account very well for charmonium and T narrow state spectroscopy and, without adjustment 
of parameters, give production cross sections which agree as well as they do with the data 
leads us to conclude that they merit further study. There are two directions for this. Exper- 
imentally, it would be interesting to compare calculated values for exclusive cross sections 
with measured values. From theoretical point of view, it may be worthwhile to evaluate 
taking predicted configuration mixing into account the allowed and forbidden Ml transition 
rates in charmonium and T and compare those with measured values. 
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